Metal-insulator transition in half-filled two-orbital Hubbard model on triangular 

lattice 
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We have investigated the half-filled two-orbital Hubbard model on a triangular lattice by means of 
the dynamical mean- field theory (DMFT). The local squared moments of charge, spin and orbital, 
and the optical conductivity clearly show that the metal-insulator transition (MIT) occurs at Uc, 
Uc=18.2, 16.8, 6.12 and 5.85 for the Hund's coupling J=0, O.OlU, U/4 and U/3, respectively. The 
distinct continuities of the double occupation of electrons, the local squared moments and the local 
susceptibility of charge, spin and orbital suggest that for J > 0, the MIT is first-order; however at 
J=0, the MIT is second-order. We attribute the first-order nature of the MIT to the symmetry 
lowering of the systems with finite Hund's coupling. 

PACS numbers: 71.30.-|-h, 64.60.-i, 78.20.-e, 71.27.+a 



I. INTRODUCTION 

Metal-Insulator transitions (MIT) and related proper- 
ties in correlated electron systems have been the central 
topics in condensed matter physics for several decadesi. 
The MIT can be easily realized by the variation of the 
external fields, doping concentration, pressure and tem- 
perature in many typical transition-metal oxides. The 
simplest and effective model to describe the low-energy 
physics of these strongly correlated transition-metal ox- 
ides is the single-orbital Hubbard model, including the 
competition between the kinetic energy and the local 
Coulomb interaction. Such a competition may result 
in many complicated and novel phenomena, such as the 
high temperature superconductivity in low-dimensional 
cuprates. Theoretically, great progress has been achieved 
in understanding the essence of the MIT in the single- 
orbital Hubbard model, mainly due to the development 
of the dynamical mean-field theory (DMFT)^ in the 
past decade. The DMFT allows us to accurately treat 
the Hubbard subbands in time axis and to obtain the 
quasiparticle peaks with the three-peak structure, which 
makes this approach an advance over the density func- 
tional theory and the Hatree-Fock approximation"^. With 
the help of DMFT, we have gotten a deep insight to many 
properties of the single-orbital Hubbard model, e.g. the 
MIT, the optical conductivity and absorption, transport, 
and so on^. Among these properties related to the MIT, 
the order of the MIT in the Hubbard model is essential. 
Bulla et ali^ demonstrated that in the single-orbital 
Hubbard model on Bethe lattice, the MIT is first-order 
for < T < Tc; whilst, it is second-order for T > Tc. In 
the two-dimensional Hubbard model, Onoda and Imadcfc^ 
also found that the MIT is first-order in finite T < Tmit 
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by means of the correlator projection approach with the 
DMFT. These suggest that in low T, the MIT in the 
single-orbital Hubbard model is first-order. 

Since the realistic transition-metal oxides, such as 
manganites, vanadates, titanates and nickelates, usu- 
ally have multiple degenerate orbitals'', the multi-orbital 
Hubbard model is more appropriate to describe the low- 
energy process than the singe-orbital Hubbard model. At 
the same time, the multi-orbital Hubbard model may 
exhibit more complicated and richer phenomena than 
the single-orbital Hubbard model. Besides the conven- 
tional localization-delocalization transition of electrons, 
there may exist many different orbital ordered phases. 
For example, in a two-orbital system, one orbital may 
be completely empty and another is fulfilled, forming 
the ferro/antiferro-orbital ordered phase; or, one orbital 
is fulfilled and insulating, another is partially occupied 
and metallic, forming the so-called orbital selective Mott 
phase (OSMP)^i^iiaiiiii^. In these situations, the orbital 
degree of freedom plays an important role in the phase 
diagram and the groundstate properties. 

More recently, a number of researches have been con- 
centrated on the nature of MIT and other properties of 
the two-orbital Hubbard modelsi^ii^ii^ii^. However, even 
on the Bethe or the hyper-cubic lattices, the nature of 
the MIT in the two-orbital Hubbard models has been 
controversial, although intensive theoretical efforts have 
already been done — i^'^'^^d^ . In the two-orbital systems, 
Inaba et al}^ and Biinemann et alr^ found that the Mott 
transition is discontinuous for any finite J > and con- 
tinuous only for J = within a generalized Gutzwiller 
approximation. However, utilizing the DMFT with the 
numerical renormalization group, Pruschke and Bulla— 
claimed that the Mott transition is second-order for J > 
U /4. They found that the variation of the local squared 
moment of spin near the transition is too small to judge 
the order of MIT for large J. By making use of the DMFT 
with self-energy functional approach, Inaba and Kog 
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believed that the nature of the Mott transition is first- 
order in all the parameter region for finite J, though they 
found that the jump of quasiparticle weight is too weak 
to identify the order of the phase transition when J is 
large enough. The controversy on the order of the MIT 
suggests that it is urgent to find a more proper quantity 
to judge the occurrence and the order of the MIT when 
J is very large. 

Up to date, most of the studies have been focused on 
the Bethe or the hyper-cubic lattices. It is not known 
what the essence of the MIT is in the multi-orbital Hub- 
bard model on the frustrated lattices. When the strong 
electron-electron interactions compete with the geometri- 
cal frustration effects, a number of unconventional phases 
and exotic properties emerge as the result of the compe- 
tition, such as the MIT and the antiferromagnetism in 
the organic compounds k-(BEDT-TTF)2X with X as an 
anio n^^i^^ , etc. Recent development in material fabrica- 
tion shows that more and more transition-metal oxides 
exhibit strong electronic correlation on two-dimensional 
triangular lattices and the multiple orbital character, 
such as NaNi02''^ and AgNi02— , etc. These appeal for 
the study on the multi-orbital Hubbard model on the 
triangular lattice. 

In this paper we focus on the MIT physics of the two- 
orbital Hubbard model on a triangular lattice by means 
of the exact-diagonalization DMFT. We adopt the lo- 
cal squared moment of charge, together with the local 
squared moments of spin and orbital, to measure the 
occurrence of MIT in the two-orbital Hubbard model, 
and find that we can well judge the occurrence of the 
MIT when the Hund's coupling J is very large. We def- 
initely show that the MIT at large J is first-order. The 
variation of the optical conductivity of the two-orbital 
Hubbard model is also consistent with the MIT with the 
increase of U. The rest of this paper is organized as fol- 
lows: in Sec. II, we describe the model Hamiltonian of the 
two-orbital system and briefly explain the framework of 
the exact-diagonalization DMFT approach; in Sec. HI, we 
present the evolutions of the densities of states (DOS), 
the local squared moments of charge, spin and orbital, 
and the optical conductivity with the on-site Coulomb 
interaction; the order of the MIT in the two-orbital sys- 
tem is also discussed in Sec. HI; the last part is devoted 
to the summary. 



II. HAMILTONIAN AND METHOD 

We start from a half-filled two-orbital Hubbard model 

(1) 



in a triangular lattice, where c\^u{ciaa) is the creation 
(annihilation) operator of the electron at site i with or- 
bital a{= 1,2) and spin (7(=tii)i ^^^d riiaa is the elec- 
tron number operator, tap denotes the hopping integral 
from the (3 orbital to the a orbital, and only the nearest- 
neighbor hopping is taken into account. For clarity and 
to compare our results with the present literature, we 
assume that the intra-orbital hopping integrals are the 
same, i.e. taa =tf3f3 =t; and we neglect the inter-orbital 
hopping, though sometimes the inter-orbital components 
may play an important role^^. 

The parameters U, U', J, and J' denote the intra- 
orbital Coulomb, inter-orbital Coulomb, Hund's and the 
pair-hopping couplings. In what follows, considering the 
realistic wavefunctions of 3d-orbitals2i and the spin rota- 
tional symmetry, we adopt the relationships J — J' and 
U = U' + 2 J. Unlike the Bethe or hyper-cubic lattice, 
the particle-hole symmetry is broken at half-filling on the 
triangular lattice. At U=0, the tight-binding dispersion 
of each orbital channel is 



eUaa -2taa[cOfi{ka;) + 2 COS ("^^y ) COS 



(3) 



with the bandwidth W ^9\t\. 

Within the framework of the DMFT, the Hamiltonian 
(1) and (2) are mapped onto an effective Anderson impu- 
rity model by integrating over all the spatial degrees of 
freedom, except for the central site o. The correspond- 
ing Hamiltonian, H^ff, contains a central "atomic" or 
"impurity" part, Hatom, and an effective medium part 
Hmed, which has to be determined self-consistently. The 
two-orbital Anderson impurity Hamiltonian reads 



(4) 



and 
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where d^^a ^^^1 al^^. create the "impurity" electron and 
a bath electron, respectively; the impurity level e^a is 
usually chosen as the zero point of energy, and the hy- 
bridization parameter Vkaa of the impurity model is cal- 
culated self-consistently in DMFT. The atomic Hamilto- 
nian Hatom of the central site is the same as H' in Eq.(2), 
and Hs represents the number of the conduction band of 
Anderson impurity model. For a set of parameters (U 
and J), we can obtain the interacting Green function 
Gaa(i^n), the free Green function GQaa{i^n) and the 
self-energy of the Anderson model as follows. 
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and the lattice Green's function is 

k 

= E[7 



double occupancy < nin2 >: 



(7) 

1/3, (8) 



and the impurity Green function G^^^ ^{iui) = 
{{daa\d^aa))^^ IS given by 



G 



I a/3 



(9) 

In Eq.(9), the spectral width function Aaa{i^) — 
'^k^kaa/i^^ + /i — ^kaa), M IS the chcmical potential, 
and ijjn = (2n + l)7r//3 is the fermionic Matsubara fre- 
quency. Throughout this paper we fix the temperature 
(3 = 1/A:bT=16. We perform the iterative procedure 
repeatedly until a self-consistent solution of the lattice 
Green's function and the self-energy are found. 

Various analytical and numerical methods can be em- 
ployed to solve the effective impurity problem. In the 
following, we first make use of the exact-diagonalization 
(ED) method to treat the impurity model Eq.(4) and (5) 
at finite temperatures of /?^. Then we perform the it- 
eration on Eqs. (6-9) repeatedly until a self-consistent 
solution. Some subroutines, such as the minimizing sub- 
routine which searches for the parameters of the Ander- 
son impurity Hamiltonian and the RS subroutine which 
is used to diagonalize the Anderson impurity Hamilto- 
nian, are from the Ref. [2]. In this paper we take ris — Q 
for each spin channel. Liebsch found that when tIs > 
3, the converged results qualitatively agree with each 
other—, and the accuracy of the ED ansatz in the single- 
orbital Hubbard model is well controlled. Demchenko et 
al. had shown that in the absence of particle-hole sym- 
metry, the pole formation and the MIT are independent 
of each other on the Bethe lattice^^-. So, the quasiparticle 
weight Z is not suitable for characterizing the occurrence 
of the MIT. In this paper, we utilize the local squared 
moments of charge, spin and orbital and the correspond- 
ing local susceptibility to characterize the nature of the 
Mott transition. 

Since the MIT is associated with the localization- 
delocalization transition of the charge degree of freedom, 
we measure the MIT with the local squared moments of 
charge, < >, 



<C'> 



together with the local squared moments of spin and or- 
bital, 



<Sl> = ((nt - nyf) 



< Tt > 



{{ni 



n2 
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All of these quantities are relevant to the spin-dependent 
double occupancy < n^ni > and the orbital-dependent 
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through a few identities, such as, 

„ OF ^BF dF. 
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here F is the free energy, and fi the chemical potential. 
The local susceptibilities of the charge, spin and orbital 
are defined by 

\T[n{T)-2][n{0)-2])dT, 

Xs = f\T[n^{r) - n^{r)][nT{0) - n^m)dr, 
Jo 

Xo = / (r[ni(T)-n2(r)][ni(0)-n2(0)])dr, (13) 
Jo 

where T is the time ordered operator, nir) — 
J2aa naair), ^^.(a) (t) = J^aia) naair), and T IS an imag- 
inary time. 

With the knowledge of the single-particle energy spec- 
trum, the optical conductivity axxi^-^) can be calculated 
in the local approximation^SiSliS^. In terms of the single- 
particle spectral weight A(k, w), a^xi^) is 



f{e)-f{e + cu) 



1 cIf 

kcr '-^ 

where e is the electron charge, $7 is the volume of the 
lattice, and e^Tr/fi is the unit of the conductivity. The 
disappearance of the Drude peak could indicate the oc- 
currence of MIT, and the optical conductivity peaks pro- 
vide many information of the charge excitation between 
subbands of the systems. 



III. RESULTS AND DISCUSSION 

As reported in the literature, the Hund's coupling 
plays an important role in controlling the Mott transi- 
tion, and the nature of the Mott transition in the two- 
orbital system on the symmetric Bethe lattice has been 
controversiaU^ii^. In the present study, we investigate 
the properties of the Mott transition on the asymmetric 
triangular lattice so as to resolve the controversial results. 

Remarkably different from the Bethe lattice, the den- 
sity of states (DOS) of quasiparticles on the triangular 
lattice is asymmetric, as shown in Fig. 1, where we also 
present the evolution of the DOS with the increase of 
the intra-orbital Coulomb interaction. For J — U/A, \t 



FIG. 1: Evolution of density of states (DOS) q{uj) with intra- 
orbital Coulomb interaction U in two-orbital Hubbard model 
on triangular lattice. From top to bottom, U=3, 5, 5.5, 6 and 



6.15; J 



(7/4, and /3 



=16.0 




FIG. 2: (Color online) Dependence of optical conductivity on 
frequency in two-orbital Hubbard model on triangular lattice. 
The intra-orbital interactions, U=3, 4, 7, and 8; J = )7/4, the 
other parameters are the same as Fig. 1. 



is clearly seen that the Mott transition has already oc- 
curred at [/ i2 6.15. A detailed numerical calculation 
shows that the critical value of the MIT is Uc=6.12. For 
other finite J, the dependence of DOS on the Coulomb 
interaction strength U exhibits similar tendency. With 
the increase of the Hund's coupling J, the critical points 
of the Mott transition occur at 18.2, 16.8, 6.12 and 5.85 
for J=0.0, O.OlU, C//4, and U/3, respectively. The ten- 
dency of Uc substantially decreasing with the increase of 
the Hund's coupling on the present 2-dimensional trian- 
gular lattice is consistent with the previous results on the 
Bethe lattice. It is very interesting that for various J, the 
critical value, Uc, of the MIT on the triangular lattice is 
about twice larger than that on the Bethe lattice. This 
may arise from two facts: one is from the spin frustration 
and fluctuation effect on the triangular lattice; another 
is from that the orbital fluctuations in the two-orbital 
system enhance the metallic character, leading to a large 
critical value, Uc. 

The optical conductivity also exhibits signatures of the 
MIT. It is interesting that how the optical conductivity 
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FIG. 3: (Color online) (a). Local squared moment of charge as 
a function of interaction strength U/W for J=0.0, O.OlU, U/4 
and U/3. (b). Dependence of double occupations < n|n| > 
(dashed line) and < nin2 > (dot-dashed line) on U for J — 0.0 
and J = 0.0117 (circle). 



evolves with the Coulomb interactions in the two-orbital 
Hubbard model on the triangular lattice. Compared with 
that of the single-orbital Hubbard model, the optical con- 
ductivity of the two-band Hubbard model is more com- 
plicated and exhibits multi-peak structure, as seen in Fig. 
2. When the Coulomb interaction U is smaller than the 
critical value Uc, the Drude peak and the charge ex- 
citation peaks exist at the same time. The multi-peak 
charge excitation structure in the present system signifi- 
cantly differs from the single-peak structure of the single- 
orbital Hubbard modeP^. The peaks a.t oj — 3.0 ~ 3.5 
and Lu = 6.0 ~ 6.5 come from the excitation between 
different Hubbard subbands below and above the Fermi 
surface. With the increasing of the Coulomb interaction, 
the intervals of these Hubbard subbands become larger 
and larger, and the charge-excitation peaks move to the 
high frequency, as seen in Fig. 2. Since the bandwidths of 
the two orbitals are identical, no orbital selective Mott 
transition is observed. As U >3, we observe a small 
low-energy mid-peak at w ~ 1.0. Such a mid-peak may 
contribute from the quasiparticle peaks near the Fermi 
level, as seen in the DOS near Ei? in Fig. 1. The exci- 
tation between the renormalizcd quasiparticle peaks and 
the Hubband subbands close to gives rise to this small 
mid-peak. When the Coulomb interaction U is greater 
than the critical interaction Uc, the Drude peak and the 
small mid-peak disappear. Subsequently, the system en- 
ters an insulating phase, as shown in Fig. 2. The insulat- 
ing gap becomes more and more wide with the increase 
of the Coulomb interaction. Up to date, the optical con- 
ductivity experiment on the compounds with two orbitals 
and triangular lattice is not available, we anticipate the 
corresponding experimental results can be done in near 
future. 

As known from the earlier literature, the Hund's cou- 
pling J plays a key role in controlling the nature of the 
Mott transition in the two-orbital Hubbard model on the 
Bethe and the hyper-cubic lattices. On the present trian- 




FIG. 4: (Color online) Local orbital squared moment vs ra- 
tio of intra-orbital Coulomb interaction U over bandwidth 
W for different J. The green dot-dashed line is for the 
SU{2)®SU{2) system with J'=0, J=0.1 and U=U'-hJ. The 
local spin squared moment vs U /W is shown for J = 0.0 and 
J = 0.01(7 in the inset. 



gular lattice, we investigate the local squared moment of 
charge and the double occupations < n^ni > (the dashed 
lines) and < nin2 > (the dot-dashed lines) on the intra- 
orbital Coulomb interaction at different Hund's coupling 
J, as shown in Fig. 3. At J/U = 0, as seen in Fig. 
3a and Fig. 3b, there is no singular jump in < >, 
< n-fn^ > and < nin2 > in this triangular frustrated 
system, which implies that the Mott transition in the 
two-orbital system with equal bandwidths is the second 
order. On the other hand, when the Hund's coupling J 
is introduced, there are discontinuous jumps in < > 
and the double occupancies, as shown in Fig. 3a and Fig. 
3b, showing that the Mott transition in the two-orbital 
Hubbard model is first-order. This result is consistent 
with what obtained via the Gutzwiller methodic on the 
infinite-dimensional hyper-cubic lattice. On the Bethe 
lattice, the jump of < 5^ > is obscure and hard to dis- 
tinguish whether the MIT is first- or second-orde r^^i^^ . 
While on the present triangular lattice, the jumps of the 
local squared moments of charge < > is very obvious 
for all J ^ 0. Therefore the MIT in the half-filled two- 
orbital Hubbard model on triangular lattice is first-order 
for all finite J situations. These results also agree with 
those of Inaba et al. by the DMFT with self-energy func- 
tional approachi?. However, at J = U/A, our result is 
contrast to that of Pruschke et al. which is obtained by 
the DMFT with the numerical renormalization group^^. 
This shows that < > is a proper measure to the MIT 
in the large J situation. 

To understand the nature of the Mott transition more 
clearly, we also calculate the local orbital and spin 
squared moments of < > and < S"^ >, as shown 
in Fig. 4. In the metallic limit of U=0, the local squared 
moments < >~< >= 1; and in the insulating 
and strongly correlated regime, < >=< S"^ >= 4/3 
for J = 0, and < >=0 and < >= 8/3 for all fi- 
nite J, which are in agreement with the linearized DMFT 



FIG. 5: (Color online) Local orbital susceptibility vs ratio of 
Coulomb interaction U over bandwidth W for J = 0.0, O.OIC/, 
C//4 and C//3. The inset shows local charge susceptibility in 
the case of J = 0.0 and J — 0.0117, respectively. 



results-^^. As seen in Fig. 4, the local squared moments 
of spin and orbital are continuous at J=0, showing that 
the MIT is the second order, in agreement with the re- 
sult from the local square moment of charge. Further, 
as seen in Fig. 4, for various finite J with J=0.01U, 
U/4, and U/3, the discontinuous jumps of the local 
squared moments of orbital and spin also demonstrate 
that the MIT is first-order, consistent with the preceding 
results. Therefore, combining the local squared moment 
of charge, < >, and those of spin and orbital, < > 
and < > , one can measure the order of the MIT over 
all of the Hund's coupling J. 

Consistent with the behaviors of the local squared mo- 
ments, the divergence of the local orbital susceptibility 
near Uc in Fig. 5, together with the discontinuous jump 
of the local charge susceptibility in the inset in Fig. 5, 
clearly shows that the MIT on the triangular lattice is 
the second order at J = 0. Since the MIT in the present 
system with finite J is the first order, the local orbital 
and charge susceptibilities also exhibit discontinuities. It 
is worthy of noticing that due to the frustration effect on 
the present triangular lattice, the local orbital suscepti- 
bility in the system with J = 0.01t7 is suppressed near 
the MIT critical point; meanwhile, such suppression is 
observed only for J = 0.03U on the Bethe latticei^. Simi- 
lar behavior is also observed in the local squared moment 
of orbital. 

Compared with the single-orbital Hubbard model on 
the triangular lattice, the critical value Uc of the two or- 
bital model much larger than that in the single-orbital 
model^i is mainly due to the orbital fiuctuations. The 
physical origin of the different order in the MIT systems 
with finite J and J=0 is still a puzzle. Biinemann et al. 
attributed it to the presence of multiple atomic energy 
scales in the two-orbital Hubbard model. This argument 
may be not true since there does exist more than one 
atomic energy scale except U in the single-orbital Hub- 
bard model. To resolve this puzzle, we suggest that the 
order of MIT in the strongly correlated Hubbard model 
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may depend on the symmetry of the systems. At J=0, 
the spin-orbital couphng system is SU(4) symmetric; on 
the other hand, the rotational symmetry of the orbital 
space is usually broken for finite J. Even if the rotational 
symmetry of the orbitals exists, i.e. U — U' + J and the 
inter-orbital Hund's coupling J' = 0, the symmetry of 
the system is SU(2)§5SU(2). However, we find that the 
phase transition in such a system is still the first order, 
as seen in the green curve (dot-dashed hue) in Fig. 4. 
We also notice that in the two-orbital Hubbard model 
with the same bandwidths, the OSMT is excluded. It 
is interesting to ask what the order of the OSMT is in 
the two-orbital triangular Hubbard models with different 
bandwidths, which deserves further study. 

One notes that in the two-dimensional triangular spin 
systems, the geometric frustration is considerable in the 
strong correlation regime, so the spatial correlations and 
fluctuations of spins may be important. In this situation, 
the approximation and precision of the present single- 
site DMFT approach should be carefully justified. For- 
tunately, when we constrain the discussion in the para- 
magnetic phases, the precision of such an approximation 
is well controlled. This has been demonstrated for the 
single-orbital Hubbard model in the triangular lattice by 
several authora^ii^S. Aryanpour et aL^ and Merino et 
al^ have shown that the results of the single-orbital 
Hubbard model obtained by the single-site DMFT ap- 
proach are consistent with those by other methods. And 
the transport properties of the 2-dimensional triangular 
Hubbard model within the single-site DMFT agrees with 
the experimental results of the organic compound'^-'^ . On 
the contrary, such a method is failed when it is applied for 
the two-dimensional square lattice. This arises from the 
fact that in the 2-dimensional triangular lattice, the spa- 
tial antiferromagnetic correlation is greatly suppressed by 
the geometric frustration, as pointed out by Aryanpour 
et al^ and Merino et alJ^ . Another reason is that the 
coordinate number of the triangular lattice is consider- 
ably larger than that of the square lattice. 



On the other hand, it is highly desirable to extend 
the present single-site DMFT approach to the cluster or 
cellular DMFT approach so as to well incorporate the 
spatial fluctuation and the intersite correlation, as devel- 
oped by many authors for the single-orbital models in re- 
cent years22i2^iMi2^i^. However, such an extension to the 
multi-orbital model meets difficulty since it goes beyond 
the ability of the high-performance computing resources 
available. And we anticipate that the cluster extension 
will not qualitatively alternate our conclusions. 



IV. CONCLUSIONS 

By using the exact-diagonalization DMFT approach, 
we have demonstrated that the Hund's coupling J leads 
to a first-order metal-insulator transition in the two- 
orbital Hubbard model with the degenerate bandwidths 
in the triangular lattice. The discontinuities of the local 
squared moments of the charge, spin and orbital show 
that the first-order metal-insulator transition occurs not 
only in the small J region, but also in the large J re- 
gion. Such distinct behaviors of the systems with finite 
J and J=0 are attributed to the lowering of the sym- 
metry of the systems. The multi-peak structure in the 
optical conductivity of the two-orbital Hubbard model 
arises from the charge excitation among more than two 
Hubbard subbands. 
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